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O , Introduction 
(D . 

Q i The purpose of this paper is to introduce and study a g-analogue of the holo- 

OO I nomic system of differential equations associated to the Belavin's classical r-matrix 
^ ■ (elliptic r-matrix equations), or, equivalently, to define an elliptic deformation of the 

>- ■ quantum Knizhnik-Zamolodchikov equations invented by Frenkel and Reshetikhin 
[FR]. In [E], it was shown that solutions of the elliptic r-matrix equations admit 
^ I a representation as traces of products of intertwining operators between certain 
^ I modules over the Lie algebra sl^. In this paper, we generalize this construction to 

I the case of the quantum algebra Uq{5lN)- 
^ ' The main object of study in the paper is a family of meromorphic matrix func- 

^ . tions of n complex variables zi,...,Zn and three additional parameters p,q,s - 
^ . (modified) traces of products of intertwiners between C/g(slAf) -modules. They are 

D I a new class of transcendental functions which can be degenerated into many inter- 
[ esting special functions - hypergeometric and g-hypergeometric functions, elliptic 
^ ^ ' and modular functions, transcendental functions of an elliptic curve, vector- valued 
^ ■ modular forms, solutions of the Bethe ansatz equations etc. 

The main result of the paper (Theorem 4.3) states that these functions satisfy 
two holonomic systems of difference equations, (4.18) and (4.19) - the first one has 
shift parameter p and elliptic modulus s, and the second one has shift parameter s 
and elliptic modulus p. 

In Section 1, we recall some properties quantum affine algebras to be used in 
subsequent sections. Our exposition follows [FR]. At the end of the section, we 
introduce the quantum analogue of the Sugawara construction for the element Lq 
of the Virasoro algebra. A similar construction is described in [ITIJMN] . 

In Section 2, we define intertwiners $(^) : Mx^k -^f./c ® 14, where Mx^k is a 
Verma module and is an evaluation representation of the quantum affine algebra 
{z G C*). Following [FR], we classify such intertwiners and prove a difference 
equation for them (due to Frenkel and Reshetikhin). Our proof uses the quantum 
Sugawara construction and is less technical than the original proof in [FR] , although 
it relies on essentially the same ideas. At the end of the Section, we reproduce the 
proof of the quantum Knizhnik-Zamolodchikov equations for correlation functions 
given in [FR]. 

In Section 3, we introduce the outer automorphism (3 of Uq{sl]\f) corresponding 
to the rotation of the Dynkin diagram (which is an A^-gon), and define the corre- 
sponding operator B acting on Verma modules. We define a class of meromorphic 
functions of n complex variables - traces of products of intertwiners twisted by the 
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equations which is a q'-deformation of the eUiptic r-matrix system. This system 
is equivalent to a system of difference equations with eUiptic coefficients, and the 
shift parameter p in this system is independent of the (multiphcative) period s of 
the elliptic function occuring in its coefficients. Systems of difference equations of 
this type were considered in [JMN], where they arose as systems satisfied by cor- 
relation functions of the 8-vertex model in statistical mechanics. It would be very 
interesting and useful to clarify the connection between the difference equations of 
[JMN] and the difference equations of this paper, but at the moment it is not clear 
how to do it. 

In Section 4, we define the fundamental trace - a matrix-valued function whose 
columns are traces of linearly independent intertwiners. We derive connection rela- 
tions for the fundamental trace which are the g-analogue of the monodromy for the 
elliptic r-matrix equations computed in [E] , and the elliptic analogue of the connec- 
tion matrices for the quantum KZ equations (see [FR]). These connection relations 
turn out to be equivalent to another holonomic system of difference equations with 
elliptic coefficients, in which the shift parameter and the elliptic modulus inter- 
change. Thus, we obtain a pair of holonomic systems of difference equations with 
elliptic coefficients such that each of them is the monodromy (= the set connection 
relations) for the other system. We call such a pair a double difference system. 

In Section 5 we describe some general properties of double difference systems with 
A^-dimensional matrix coefficients. We start with conventional holonomic systems 
and show that a nondegenerate system with rational coefficients is always solvable 
in meromorphic functions, whereas for systems with elliptic coefficients in more 
than one variable this is not true. However, if a system with elliptic coefficients 
is solvable, and one fixes a matrix solution of it, then one can construct a dual 
system satisfied by this solution, also with elliptic coefficients, which is a system of 
connection relations for the previous system, and in which the shift parameter and 
the elliptic modulus interchange. In this way we obtain consistent double difference 
systems with elliptic coefficients. We give a complete classification of such systems 
in one variable and one dimension, and compute their solutions. This result can 
be extended to many variables, but in more than one dimensions the classification 
is hardly possible. Even the problem of finding a dual system to a given difference 
system with elliptic coefficients (so that they make a consistent double difference 
system together) is a difficult trancendental problem (which is no surprise since it is 
a generalization of the notoriously difficult problems of finding connections matrices 
for difference systems with rational coefficients, and of computing monodromy of 
differential equations). In fact, the example given in this paper appears to be the 
first nontrivial explicit example of a consistent double difference system with elliptic 
coefficients. 

In the appendix, we briefly discuss some limiting cases of the difference equations 
deduced in the preceding sections and their connections to qKZ equations, elliptic 
KZ equations, Smirnov's equations, and Bethe ansatz equations. 

Acknowledgements I would like to thank my advisor Professor Igor Frenkel for 
formulating the problem and stimulating my work, and also I.Cherednik, J. Ding, 
M.Jimbo, M.Kashiwara, D.Kazhdan, T.Miwa, Ya. Soibelman, L.Takhtajan, and 
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1. Quantum afRne algebras. 

Let 9 be a finite dimensional simple Lie algebra over C of rank r. Denote by 
<, > the standard invariant form on q with respect to which the longest root has 
squared length 2. 

Let [) denote a Cartan subalgebra of g. The form <, > defines a natural identi- 
fication [)*^f): Ai-^/iAforAGf)*. We will use the notation <, > for the inner 
product in both i) and ()*. For the sake of brevity we will often write A instead of 
hx. 

Let CKj, 1 < z < r be the simple positive roots of Q. Let 9 be the highest 
root of - the positive root such that 9 + ai is not a root for any i. Extend the 
Cartan subalgebra [) by adding a new element c orthogonal to f). Denote the Cartan 
subalgebra extended by c by f): = Cc. 

Introduce the convenient notation ckq = —9, Hi — ha^, 1 < i < r, Hq — —he + c. 
Let ai, = ^^,0<i,j<r. 

Let p be the half sum of positive roots of g, and let A?" = 1+ < p, ^ > be the dual 
Coxeter number of g. 

Let t be a complex number and q = e^. We assume that |g| < 1. If A is a number 
or an operator, by q"^ we will always mean e^^. 

Let Uq{g) be the quantum affine algebra corresponding to g [Drl,J]. This algebra 
with unit is generated by elements e^, /i, < z < r, q^'^, h E t), satisfying the 
standard relations 



q'^Biq-" = q'^'^^^Bi, q^fiq'" = q"^'^^^ h , < z < r, 

qHi _ q-Hi 
(^ifj - fj^^i = % : < i, j < r, 

n=l ^ ^ 9i 



n= 

1-a 



n=l ^ ^ "Ji 

(LI) 5i = g2/<«.-«.>, 0<i,j<r,i^j, 



771/ 

where ( ) is the g-binomial coefficient: by definition, 



n 



(1-2) 



The comultiplication in Uq{g) is defined as follows 



A{q^)^q^^q\ hei), 
A{ei) =ei^q^' + l^ei 
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The counit is given by 

(1.4) e(e,) = e{fi) = 0, < i < r, e(g^) = 1, /i e ^. 
The antipode is given by 

(1.5) S{ei) = -ei5-^^ S{fi) = -q^'^fi, < i < r, 5(g'^) = g-^ hei 

Equipped with aU these structures, Uq{Q) becomes a Hopf algebra. 
Let us add a new element D to Uq{g) satisfying the relations 

[D,ei] = a, [DJi] = -fi, 0<i< r, [D,h] =0,hel) 

(1.6) A{D) = D®1 + 1®D, e{D) = 0, S{D) = -D. 

The element D defines a grading on Uq{g): an element a is of degree n if [-D, a] = 
na. An element is called homogeneous if it is an eigenvector of the operator of 
commutation with D. 

We denote the Cartan subalgebra {) extended by D by 1^, and the Hopf algebra 
Uq{Q) extended by D by Uq{g). 

Let us define two classes of modules over Uq{Q): Verma modules and evaluation 
representations (our definition is slightly different from the standard one; cf. [E]). 

The definition of a Verma module of highest weight A and level k is slightly 
twisted: Ma,a;, (A e i)*ik e C), is the module freely generated by a vector v 
satisfying 

(1.7) BiV = 0, qf'v = q<^+kp,h>y^ h e i), q^v = q^^v, Dv = ^j^^^v. 

2(/c -|- 1) 

Evaluation representations are defined as follows. Let y be a representation 
of Uq{Q) which belongs to category O when restricted to the subalgebra Uq{Q) 
generated by ei,/i i > 1, q^, h e f). Consider the representation of Uq{g) in the 
space of F-valued Laurent polynomials V (g) C[z, z~'^], defined by 

(111) ( Zi 

(1.8) Dow{z)=z — - — -; a o w{z) = z'^77v{a)w{z), aeUq{Q), deg{a) — n. 

dz 

This representation is reducible. Let V{z) denote the submodule in it spanned by 
the vectors v(>Siz\ where v E V runs over all vectors of such weights A that < A, p > 
—I is divisible by N. Then V{z) is called the Laurent polynomial representation of 
the quantum affine algebra associated to V. 

If zq is a nonzero complex number, then the space of all elements of V{z) vanish- 
ing at Zq is a Uq{Q)- submodule of V{z). The corresponding quotient is denoted by 
Vzo- We will say that Vz^ is the evaluation representation (at 2:0) associated to V. In 
particular, if = 1 then V^q = V, so V itself is also an evaluation representation. 

The universal quantum i?-matrix for Uq{g) is defined as follows. Consider the 
subalgebras Uq{n^) and Uq{n~) in Uq{g) generated by {cj} and {/i}, respectively. 
It is known [Drl] that there exists a unique pairing between these two algebras, 
<, >t such that < Cj, fj >t— t~^Sij, and 

< a6, c>i= ^ < a, n >t<b,pi >t, A{c) = ^pj iS)rj,a,b e Uq{n'^), ceUq{n~), 

j j 

(1.9) 

< c, ah >t= ^<pj,a >t< Tj, b >t, A(c) = ^pj® rj, a,be Uq{n~), c e ?7g(n+) 
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Let {tti} be a basis of Uq{n'^) consisting of homogeneous vectors, and let {a*} be 
the basis of Uq{n~) dual to {ai}. Then the universal i?-matrix is given by 

(1.10) n = g'^^d+d^c+E •=! xj(Sxj J2ai® a\ 

i 

where {xj} is an orthonormal basis of f), and d = jj{D — p). 

The J?-matrix should be regarded as an infinite expression which makes sense 
as an operator in the tensor product of Verma modules. Such expressions form an 
algebra - a completion of the tensor square of the quantum affine algebra, and the 
i?-matrix is its element. 

The i?-matrix satisfies the following quasitriangularity axioms: 

nA{x) = A°p{x)n, 

(1.11) (Id® A)(:^) = :^i3:^i2, 

where A"p denotes the opposite comultiplication and 

i i 

iZi2 = (g) 6j (8) 1 = (g) 1, 

i 

i 

(1.12) TZ23 = '^l®ai®bi = l^TZ. 

i 

The matrix R also satisfies the quantum Yang-Baxter equation: 

(1.13) ^^12^13^23 = ^23^13^12 

We introduce the modified /^-matrices 

(1.14) 7^ = q-C®D-D®c^^ ^1^^^ ^ ^^D ^ X)qi(^z-D ^ ;l), C = 

Note that TZ{z) is a power series in z which includes only nonnegative powers of z. 

The matrix TZ satisfies a modified version of the quasitriangularity axioms and 
quantum Yang-Baxter equations: 

7^A(x) = Ag^(x)7^, 

(A®Id)(7^(^)) =7^f3(^)7^23(^), 
(Id (g) A)(7^(^)) = n-{^{z)ni2{z), 

(1.15) 7^l2(2;l/2;2)7^f3(2;l/2;3)7^23(2;2/2;3) = 7^23(-^2/-^3)^^3'^(-^l/2;3)7^l2(-^l/-^2), 
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Unlike TZ, the element TZ has a remarkable property: it can be projected to eval- 
uation representations. Namely, for any z & C* and an evaluation representation 
Vz one can define the elements 

L+{z) = (Id®7rvJ(7e°P) = (Id(8)7ry)(7e°P(^)) e Ugig)^EndiV), 

(1.16) Ly{z) = (Id®7rx/J(7^-^) = (ld^7Tv)in-\z-^)) e t/g(g)®End(y), 

where ® denotes a completed tensor product. These elements are called quantum 
currents. 

If we pick a basis in V labeled by a set / then the quantum currents can be 
viewed as matrices, L^, i,j e /, whose entries are Laurent polynomials of z with 
values in Uq{g). 

Further, one can define the projection of TZ into two evaluation representations 

(1.17) R^'^\zi, Z2) ^ {ttv^ ®TTv^)(n). 

This projection turns out to be a power series in 2 = z\j z^ with only positive powers 
of z present. This series converges in a neghborhood of the origin and therefore 
defines a holomorphic function with values in End(F-^ (8)F^) in the neighborhood of 
the origin, which we will write as i?^ ^ {£) or, when no confusion is possible, simply 
as R{z). One can show that this function extends to a meromorphic function in C 
which is a product of a scalar trnscendental function (p and a rational matrix- valued 
function R: R^ ^ (z) = (j)^ ^ iz)R^ ^ (z). The rational function R regarded as 
a function of y — log 2; is a trigonometric solution of the quantum Yang-Baxter 
equation. It satisfies the unitarity condition R{z)R2i{z~^) = 1. 

Let us define a product operation applicable to quantum currents. Let a = 
ai (g) a2 e Uqio) (g) End(y), 6 = 61 (g) 62 e Ugio) (g End(II^). Define a "product" of a 
and 6 by a * 6 = aibi (g 02 (g 62- It is important to distinguish this operation from 
the usual product. 

Let us now write down the commutation relations for currents. 

Proposition 1.1. [PR] The following relations between power series with values 
in End{Mx,k (g 1^^ (g V^) hold true: 

(1 ® R{-))L+, (zi) * L+, (^2) = Svw^ {L+2 {Z2) * L+, {zi)){l ® i?(-)), 

Z2 V V ^2 



(1 ® R{ — ))L-^{zi) * L-2iz2) = Svw^{L-2{z2) * L-,{zi))(l i?( — )), 

Z2 Z2 



(1.18) 



— k k 

(1 (g R{^-^))L+, (zi) * {Z2) = Svw^ {L-2 {Z2) * L+, {zi)){l g) R{^)), 

Z2 Z2 



where Sy^v^ denotes the permutation of the and factors. 

To prove this proposition, it is enough to apply the maps TTyi (g ttv^ <g ttmx k ' 
T^Mx k ® '^v^ ® '^v^ 1 ciiid TT^i (g ttmx k ® to the quantum Yang-Baxter relation 
for TZ. 

Let us now describe a quantum analogue of the Sugawara construction. Let 
m : Uq{Q) g) Uq{Q) — > Uq{Q) be the multiplication map. Consider the element 

(1.19) u = m{{S® Id) (n^'P)). 
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Proposition 1.2. [Dr2] 

(1.20) u-' = J2S-\b*)al where n-^ = Y,<®K\ 

(1.21) UXU-^ = S'^{x), X e Uq{Q); 

(1.22) A(w) = (w 'u)(:^°p:^)-i = (:^°p:^)-1(w ® u). 

Thus we have the foUowing proposition. 
Proposition 1.3. In the Verma module Mx,k 

(1.23) u = q^^. 

Proof. First of all, we have the equality q^^xq~^^ — S'^{x). This equality can be 
easily checked: it is enough to check it for the generators x = Si, fi,q^, since both 
sides of it are automorphisms of the quantum affine algebra. Hence, it follows from 
the previous proposition that uq~^^ commutes with the quantum affine algebra, 
so it is a constant. To prove that this constant is 1, it is enough to check that 
uv = q^^v, where v is the vacuum vector, which is straightforward. 

Let 

(1.24) n = (^-'^^ ® l)7^(Q2c^ ^ 1). 

Proposition 1.4. ( Quantum Sugawara construction) The following relation is sat- 
isfied in any Verma module: 

(1.25) ^2(c+i)D ^ ^ Id){n°P)) 

Proof. It follows from the definition of it that the right hand side of (1-54) is equal 
to q^'^^u. But u = q^^, so we get (1.25). 

2. Intertwining operators and difference equations 

We will be interested in C/q(0)-intertwining operators ^{z) : Mx^k Mjy^k^z'^V{z).i 
where (8) denotes the completed tensor product, and A is a complex number. It 
turns out that such operators may be nonzero if and only if A equals 2('fc+i)^''^^ 
plus an integer. The shift by an integer is unimportant, so we will assume that 

A <i^,i^> — <A,A> 

^ - 2(fe+l) 

Proposition 2.1. [FR] Operators $ are in one-to-one correspondence with vectors 
in V of weight X—v. This correspondence is defined by the action of^ at the vacuum 
level. 

Let zq be a nonzero complex number. Evaluation of the operator ^{z) at the 
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Sometimes (when no confusion is possible) we will use the notation $(2) for the 
operator $ evaluated at the point z E C*. This will give us an opportunity to 
regard the operator ^{z) as an analytic function of z. This analytic function will 
be multivalued: ^{z) = z^^'^{z), where is a single- valued function on C*. 

Let W be an evaluation representation of Uq{Q). Then the intertwining property 
for ^{z) can be written in the form 

mz) ® Id)L+ («;) = Rwv{q~''w/z)L+,i'u^)mz) (g) Id) 

(2.1) {<!>{z) ®ld)L^{w) = Rvw{z/w)-^L^iw)mz) ^Id). 

Relations (2.1) are equalities of maps Mx^k'^W — > Mi.^k'^Vz ® W. 

To prove these formulas, it is enough to combine the second and third relations 
of (1-15) with the intertwining relation ^{z)x = A{x)^{z). 

Remark. Note that relations (1.18) and (2.1) are a priori satisfied only 
formally, as equalities between power series. However, since we know that both 
sides of these equalities extend to mcromorphic functions, we can conclude that 
they are also satisfied analytically for almost all values of the parameters. 

Let us now deduce the difference equation for intertwining operators, following 
the method of Frenkel and Reshetikhin. 

Let 

(2.2) U = q^^^'u. 

Introduce the notation Q — g-2(fc-i-i)-D^ ^ _ ^-2(fc-i-i)_ assume that \p\ < 1. The 
quantum Sugawara construction implies that Q = in M\^k- Also, the operator 
Q acts in the Laurent series representation V{z) as follows: Qv{z) = v{pz). 

Since ^» is an intertwiner, we have the following relation between Laurent series 
in z: 

(2.3) Hz)Q-'^ = A{Q-^)^{z) = (Q-^ ® Q-^)$(z). 

Using the quantum Sugawara formula (1.25) and the identity (1 iS> Q~^)^{z) = 
^{p~^z), we obtain 

(2.4) ^p-^z) = {Q® l)^z)Q-^ = (U-^ l)Hz)U, 

where U is defined by (2.2). 

Introduce the notation $»(XljOi®^i) = I]j(l®^i)*'^i: ^ ^ }iomc{Mx,k, My^k<^z^V{z)). 

Lemma 2.2. 

(2.5) ® l)^z)U = (L+(g^^)-i$(^)) . Ly{p-'z). 

Proof. We have U — uq^'~"^ = S{bj)ajq^'~'^ . Therefore, since ^{z) is an inter- 
twiner, 

(2.6) ^z)U = Y,^{S{bjmz)ajq'^^ = ^(5 ® 5)(A°P(6,))$(^)a,g^^^. 
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Following Drinfeld ([Dr2]), introduce the notation {X <^Y <^ Z) o ^ = {S{Y) (g) 
S{Z))^X. This defines a right action of the tensor cube of the quantum affine 
algebra on the space Homc(MA,fc, Miy^k^z^V{z)). 

Using this notation, we can write (2.6) as follows: 

(2.7) ^z)U = (Id ® A°P)(^) o $(2) . g2C"^. 
Applying (1.11), we get 

(2.8) <^{z)U = (^12^13) o ^z) ■ = o {ni2 o $(^)) • g2Ci?. 

Let us separately consider the expression X = iZi2 o ^(z) which occurs in (2.8). 
Using the intertwining property of ^(z) and (1-11), we obtain 

(2.9) X = Y^z), Y = m3i((A ® S){TZ)), 

where 777,31 {a <S> b <S> c) = ca <S> b. 
Applying (1.11), we find 

(2.10) Y = m3i ((Id ® Id ® 5)(:^i3:^23)) = {S^ Id)(:^°P) • {u 1). 
Thus we have 

{U-^ l)^z)U = {U-^ (8) 1) [iZi3 o{S^ Id)(:^°P)(w l)^z)] q'^^^ = 

(using (1.21)) 
(9-2^^ 1) [iZi3 o (S-^ ® Id)(:^°P)$(^)] g2C?^ = 

(g-2^^ ® l)[^i3 o (^°p)-i$(^)]q2C^ = 
(g-2^^ ® 1) [(^°P)-i$(2) • (Id ® 5)(^i3)]q2^^ = 

(q-2^^ ® 1) [(^°P)-1$(z) • (S--^ ® Id)(^r3')]9'^'' = 

(g-^'^^ 1) • (7e°P)-^(l (8) g-'=^)(g2'=^ g2'=^)$(z) • {n^^{q-^-^''){l (8) g-'^^)) = 

n^,\q^Mz).n^,\p) = 

(2.11) (L+(g'=^)-i$(^)).L-(p-i^). 

The lemma together with equation (2.4) implies the following difference equation 
for ^{z): 

Theorem 2.3. [FR] The intertwining operator ^(z) satisfies the difference equa- 
tion 



(2.12) ^pz) = L+{pq^z)mz) • L-{z)-'). 

Let us now deduce the difference equations for quantum correlation functions. 

Let V^^.-.^V^ be evaluation representations of Uq{Q), and let ^^{zj) : M\.^k 
Mx-_^^k®yij be t/g(0)-intertwining operators. Then the product <^^{zi) . . . $"^(2^) 

™„i „„ J- — i\/r . i\/r x^T/l ^ ^ T/n ;f I ~ I ^ ^ 1^1^^ 
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••• >> \zn\. Consider the matrix element of this product corresponding to the 
vacuum vectors in the Verma modules: 

(2.13) zn) =< vl, ^\zi) . . . $"(^nK„ >, 

where va„ is the highest weight vector of Mx^ and is the lowest weight vector 
of Ma(, . This function takes values in the tensor product - ■ • ^F". It is called 
the quantum correlation function. 
We have 

^{zi, ...,pzj..., Zn) =< vl^,^\zi) . ..^\pzj) . . .$"(2;„K„ >= 

(2.14) 

< vl,^\z,) . ..^^-\zj-^)L+,{pq>'zj)^\zj) . L-,izj)-'^^+'izj+,) . . .$n^n)^A„ > 

Let us now drag to the left and L~ to the right, using commutation relations 
(2.1). Taking into account the relations 

{S Id){L+{z)){vl ^v)= vl q^°v, 

(2.15) L-{z)-^{vx^ ^v) = vx„ ® q^-v, 

which follow from the definition of the quantum currents, we will get the following 
result. 

Theorem 2.4. [FR] The quantum correlation functions satisfy the following sys- 
tem of linear difference equations. 

^{Zi, ...,pZj...,Zn) = ti,j_-i^ { ( " )\Vj>^ 

Zj-i Z\ 

(2.16) { — ) ...Kj+ij {^—) ^{zi,...,Zj,...,Zn). 



3. Difference equations for traces of intertwiners. 

From now on the letter q will denote the Lie algebra sZjv(C) of traceless N x N 
matrices with complex entries. The dual Coxeter number of this algebra is N, and 
the rank is N — 1. The Cartan subalgebra [) is the subalgebra of diagonal matrices. 

Let B be the N x N matrix of zeros and ones corresponding to the cyclic per- 
mutation (12.. .iV): 





( ° 


.. 


. 


1 


\ 




1 


.. 


. 







B = 





1 .. 


. 









V o" 


.. 


1 





/ 



Define an inner automorphism (3 of f): /3(a) = BaB~^, a G f). This automorphism 
has order N. We can extend it to 1) by making it act trivially on c and D. 

The automorphism /3 can be extended to an outer automorphism of Uq{Q) defined 
by the relations 

1 ^ or^ \ Ql f \ — f Q{ „h\ B(h) u ^ Z 
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where the subscripts are regarded modulo N. This outer automorphism also has 
order N and can be interpreted as the rotation of the Dynkin diagram of g (which 
is a regular A'^-gon) through the angle 2'k/N . 

The action of 13 in Uq{Q) preserves degree, hence, it preserves the polarization. 
Therefore, it transforms Verma modules into Verma modules. In other words, we 
can regard (3 as an operator B : Ma,^ — > Mp(^x),kj where by convention f3{\){h) = 
X{f3~^{h)). This operator intertwines the usual action of Uq{Q) and the action 
twisted by (3: [3{a)Bw = Baw, a G Uq{Q), w G M\^k- 

Let u = /3~^(A), and let (^^{zj) be as above (cf. section 2). 

We assume that the representations are finite-dimensional and irreducible 
when restricted to Uq^Q). It is easy to show that for any such representation V there 
exists a unique, up to a constant, operator B -.V ^ V such that Bav = (3{a)Bv 
for a G C/g(0), f e y - it follows from the fact that V twisted by (3 is isomorphic 
to V. The operator B gives rise to a well defined automorphism (3 of End(V): 
(3{E) = BEB-^. 

Let s be a complex number, < |s| < 1. Following the idea of Bernard [Ber], 
Frenkel, Reshetikhin ([FR], Remark 2.3) (see also [E],[EK]), introduce the following 
formal power series in Zi/z2, z^j sz^/zi. 

(3.2) F{zu zn\s) = TV Im,,, (^'(-^i) • . . $"(^n)Ss-^). 

It is not difficult to prove that this series defines an analytic function when \zi\ » 
\z2\ » ■ ■ ■ » \zn\ » \szi\. This function takes values in F-*^ (8) • • • (8) F". From 
now on it will be the main object of our study. 

It turns out that the n-point trace F{zi, Zn\s) defined by (3.2) satisfies a 
remarkable system of difference equations involving elliptic solutions of the quantum 
Yang-Baxter equation for sIn- Let us deduce these equations. 

According to (2.14), we have 

F{zi, ...,pzj..., zn\s) = Tr |m,,, {^\zi) . . . ^\pzj) . . . $"(^n)Ss-^) = 

(3.3) 

Tr Im,,, {^\zi) . ..<^'-Hzj-i)L+,{pq'^z,)^^{z,) . L-,{z,)-'^^+\zj+,) . ..^-{zn)Bs 
We would like to deduce a difference equation for F. 

First of all, we need to describe some properties of the automorphism p. Since 
/3 is a degree preserving automorphism, it preserves the universal i?-matrix: 

(3.4) {(3 ^ p) {TV) = n, (/3 ® p) (7^) = 7^. 

Let = iP^ (X) l)(7^). It is obvious that {P^ ® P^){T^) = ■ 
Introduce the following notation: 

L+(^)^ = (Id®7^^.J((7^^)°P) 
L^(^)^ = (Id®7^vJ((7^^)-^), 

(3.5) B^"''\z)m = i.'KvHz) ® 7rvJ(7^^). 

Let be a finite dimensional representation of Uq{Q). Let $ : Mi3(^x),k ~^ 
M\^k <^W he a,n intertwining operator. Introduce a new function F^"^ {x,y\s) with 
values in (g) EndF EndF defined by 
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Our plan is to show that F satisfies a system of difference equations, by dragging 
L"*" around the circle from right to left and L~ from left to right using the commu- 
tation relations for quantum currents and the property of trace: Tr(afe) = Tr(6a). 
Before we do so, we need a few identities. 

First, we have the following generalization of the third equation in (1.18): 
(3.7) 

L+,{zi)iR32{q''zi/z2)j+i{Ly,{z2)i) ^ = {Ly,{z2)i) ^R32{q Zi/ Z2)j+lL+,{zi)l 

which implies that 
(3.8) 

iL-,{z2)i)-'*L+,{z^)i = Rli{q''z^/z2)j+iR'^2{Q^^ 

where R^'^ implies that the first component of R is applied from the right and the 
second one from the left, and R^'^ implies that the first component of R is applied 
from the left, and the second one from the right. These relations are equalities 
between elements of the product Uq{g) (X> Endy Endy. The notation L{z)2 and 
L{z)s implies that the second component of L{z) operates in the second and third 
factor of this tensor product, respectively. To prove these relations, it is enough to 
apply the automorphism 1 (8) to the quantum Yang-Baxter relation for 
the i?-matrix, and then project the obtained relation to the corresponding product 
of representations of the quantum affine algebra. 
Next, we have the identities 

(3.9) Lyizrh-"" = s-^'Lyis-h)-', s-^L+{z) = L+{sz)s-'', 
and 

(3.10) L-{z)-'B = Bp-\L-{z)-'), BL+{z) = P{L+{z))B. 
Finally, we have the following /3- twisted versions of relations (2.1): 

^iL+{x)i = R3i{q-''x)iL+{x)i^i, 

(3.11) {L-{y)i)-'R^2{y-')I^l = ^i{L-{y)i)-\ 

These relations are equalities of elements of Hom(M;\^fc, Mi^^k^W) ^EndV <SiE'n.dV , 
and the meaning of the subscripts 1,2, and 3 is as above. 

Now we are in a position to compute F. Combining relations (3.6)-(3.8) with 
(2.1), we get 

F'+''{sx,y\s) = i2S2(^)j+/+i^32( — )7i/+i^3i(p^x),+iF^^(x,|/|.), 

y y 

(3.12) 

y iy iy 

Let Q^^ = B(B~^F^'^B~-^Bl. Then we have 

g'+^\sxMs) = RU^-^)Rz2{^)~^R^i{psx)B^'g"{xMs)B3. 

(3.13) 
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Let 

y y 

(3.14) T^ix^y) = S^<2(^)42(^)i?5^(^)-^(S2)-\ 

and let 
(3.15) 

Po{x,y) = Toix,y)To{sx,y) . . .To{s''-^x,y), P^{x,y) = T^{x,y)T^{x, s'^y) . . .Ti(x, 5"^+^), 

where Tq, Ti, Pq, Pi e End{W ® EndV EndF). Let be denoted by G. Then 
we have the foUowing equation: 

oo oo 

(3.16) g{x,y\s) = n ^o(x,^^^y) n ^^1(0,5-^^^)^^(0,001^) 

M=0 M=0 

Let us now consider the expression ^(0, oo|s). We have 

^(0, oo|s) = 

(3.17) 

Apart from that, we have the identity 

TV Im,,, (/i$iSs-^) + /iTr Im,,, ($iSs-^) = 

(3.18) Tr |m,,, ($i/iSs-^) = IV |m,,, (/^-H/i)*!^^-'') , e I), 

which implies that 

(3.19) IV Im,,, (/i^iBs-^) = (/3-1 - l)-i(/i)IV Im,,, (^iSs-^). 
Hence, if is any analytic function then 

(3.20) IV |m,,, (0(/i)$ii?5-^) = (/>((/3-^ - 1)-'(M)IV |m,,, (^iS^-''). 

For brevity introduce the notation x — (/?~^ — Then expression (3.17) can 

be rewritten in the form 

^(0, oo|s) = 

(3.21) 



l®{kp®l-l®kp)+'£^{x{xi)®l^Xi-{l+x){xi)mi 



IV Im. . ($i5s-°) ® Id ® Id 
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Let 

P{x,y) = 

(3.22) 

oo oo 

Y[ Po{x,s^^y) Y\ Pi(0, s~'^-'^y)5^®('^'°®^~-^®'^''^+^i(^(^*)®^®^*~^^+^)(^*^®^*®^) 

M=0 M=0 

In our situation W = Vi ® V2 ■ ■ ■ ® Vn, V = Vi, and EndVi naturally acts on 
W. Therefore, it makes sense to consider the EndVF-valued function 

(3.23) X,{z^,...,Zn\s) = rm2i{P^'^{z,z)), 

where P^^{x,y) e EndW ® EndF^ ® EndF^ is defined by (3.22) with V = Vi, and 
by definition 771321(0 ® 6 ® c) = cba. 

Now it remains to observe that equations (2.12), (3.16), (3.21), and (3.22) imply 
the following difference equation for the function F: 

Theorem 3.1. The function F{zi, Zn\s) satisfies the difference equations 

(3.24) F{zi, ...,pzi, Zn\s) = Xi{zi, Zn\s)F{zi, Zi, Zn\s) 

Note that the coefiicients of these difference equations are meromorphic functions 
which we have explicitly represented as (contracted) infinite products of trigono- 
metric i?-matrices. In the next section we will show that they can be expressed in 
terms of elliptic functions. 

4. Monodromy equations 

The quasi-classical limit (i.e. the limit — > 1) of equations (3.24) is the system 
of elliptic KZ equations described in [E] . For this system one can define the notion 
of monodromy which turns out to be expressed by products of i?-matrices. It is a 
natural question what is the quantum analogue of monodromy. In this section we 
will give an answer to this question. The answer is that the role of monodromy is 
played by another system of difference equations with elliptic coefficients which are 
products of /^-matrices depending on spectral parameters. 

Let us first describe how to interchange the order of intertwining operators. 

Let ^'^'-^''^(z) : M\^k My^h®z^V {z) be the intertwining operator such that 
< f *, $™''^''^(2)fA >= w, w e V^~'^. Suppose that ^1,2:2 are nonzero complex 
numbers, and we have a product $«'i.^i.^o(2;i)$^2,A2,Ai(^^2^ . Mx2,k Mxo,k<^V'^i^ 
where and V"^ are finite dimensional representations of Uq^Q). The question 
is: can this product be expressed in terms of products of the form $(22)$(-2i)? Of 
course, we can only talk about such an expression after analytic continuation, since 
the former is defined for \zi\ » \z2\1 and the latter for \zi\ « \z2\. However, if 
we apply analytic continuation, the answer to the question is positive, and given 
by the following theorem. 

Theorem 4.1. (see [FR]J Let Xii, he a basis of {V^)'^~'^° , and let yii, he a basis of 
(y2)A2-i^. r/ien 

j^V^V^ ^^yl^Xr-x^ ,Ai,Ao ,A2,Ai ^^^^ ^ 

(4.1) A J2 ^,^:.',\v(-)''''''^^^*"''"'^°(-22)^"-'^^'"(-2i), 
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where A is the analytic continuation, E^'t^ is a matrix called the matrix of exchange 
coefficients, and P is the permutation: ^V"^ ^ V"^ . 

Clearly, the matrix E^'^^{z) (we drop the superscripts V^,V^ when no confusion 
is possible) represents a linear operator (V^ <SiV^)^~^ — > {V^ <SiV^)^~^^. Therefore, 
if we define 

(4.2) E^{z) = ®^E'''^'{z), 

then E^(z) will correspond to an operator: (E) ^ V'^ ^ . 
Let us introduce some convenient notation. Define the operators 

Ej{z)v\...,v^ : y ^ ® ■ ■ ■ ® ® V^+^ • • ■ ® ^ 1/^ ® ■ ■ ■ ® V^+^ ^V^^-'-^V 

as follows: if Vi e V^, 1 <i < n, and hvi = Xi{h)vi, Xi G f)*, ^ G f), then 
(4.3) 

Ej{z)v-^,...y^{vi®- ■ ■(^Vj(S)Vj+i®- ■ -^Vn) = fi®- • ■(S)E^^+'-{z)^'^' {vj(^Vj+i)(^- ■ -(^v, 
where Xj are defined by 

(4.4) A,- = iP- 1)-' (j2 ^0 + E ^ -^^ ^ ^- 

Vi=i / i=i 

Let tj denote the elementary transposition (jj+1) in the symmetric group S^- 
For every a G Sn set Ej{zY — Ej{z)y<y{i)^ y,T(n). This operator maps V'^'^^^ ® 

■ ■ to . (we make a convention that for two permutations 

<7i, CF2 o'iO'2{j) — Ci(c2(i)), 1 < i < i-e. the factors in a product of permutations 

are applied from right to left). 

Let us describe some properties of the matrices of exchange coefficients. First of 

all, Ej{zY is not a single- valued function of z, so one should consider the function 

B,{CY = E,{e^Y. 

Proposition. [FR] The functions Bj{Q have the following properties: 

(i) double periodicity: 

(la): B,{C + \ogpY = B,{CY; 

(4.5) (ih): BjiC + 27rV^Y = LBj{CYL', (logp = -2t{k + 1)) 

where L, L' are constant matrices; 

(a) the quantum braid (Yang-Baxter) relation: 
(4.6) 

^,(Cl-C2)"i3, + i(Cl-C3)*^i3,(C2-C3)*^ + ^*^- =^, + l(C2-C3)"i3,(Cl-C3)*^ + ^i3, + i(Cl-C2)^ 

(Hi) unitarity: 
(4.7) B,{CY%{-CY^ =1. 

Remarks. 1. Properties (ib), (ii) and (iii) follow directly from the definition 
of Bj. Property (ia) follows from the fact that matrix elements of products of 
intertwiners satisfy the quantum KZ equations. 

2. Statement (i) implies that the matrix elements of Bj must express in terms of 
elliptic functions. Frenkel and Reshetikhin showed that in some special cases these 
functions are the elliptic solutions of the quantum Yang-Baxter equation which 
occur in statistical mechanics. 

From now on we will be assuming that the matrices of exchange coefficients are 
known. 
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Definition. The fundamental trace is a function T of zi, Zn,p, s, q with values 
in the space End{V^ (8) • • • (8) V^) defined by the property: if v = vi <S> V2 <S> • • • <S> Vn is 
a vector in <S> ■ ■ ■ <^ , Vi e V^, 1 < i < n, and hvi = Xi{h)vi, Xi & ^* , h & i), 
then 

(4.8) J^v= IV Im;,,^ ($''^'^^'^°(^i)...$''"'^"'^"-n-2n)5s-^), 

where Xj are defined by (4-4) ■ 

It follows from Section 3 that the fundamental trace satisfies p-diffcrence equa- 
tions (3.24). Below we will show that it also satisfies another system of s-difference 
equations. 

Indeed, let us carry the operator ^{zi) from left to right. We have to interchange 
it with $(^i+i), . . . , ^{zn), then with B and (this will change ^{zi) to (1 (8) 
B)^{s~^Zi)), then to shift it to the left side of the product, using the property 
Tr(a6) = Tr(6a), and then interchange it with $(-2i), • . . , ^{zi-i). This procedure 
allows us to deduce a difference equation for J^. Before we write it down, let us 
introduce some notation. 

Let ajm = tm-i ■ ■ ■ tj+itj, j <m<n, ajm = tm ■ ■ ■ tj-2tj-i, 1 <m < j. 

Theorem 4.2. The fundamental trace T satisfies the relations: 



i?]^lj^'(^^) ^T{z\, ...,Zj, Zj+i...,Zn\s) = F{zi, ...,Zj+i,Zj, ...,Zn\s)Ej{—^) 
Zj Zj+i 




Id 



where Bj denotes the action of B in . 

Now consider the system of difference equations 

G(zi, S ^Zi, Zn) = Mj(zi, Zn)G(zi, Zi, Zn), 

(4.10) 

M,(.„ Zn) = iqjT'i^) . . .<^^(^) v<^^(7^)-^ • • 

J — 1 1 j j 

This is one of the versions of the quantum KZ equations. Let Go{zi, Zn\s) denote 
the fundamental solution of this system of equations - the solution with values in 

log Mj 

End(yi ® ■ • • y") such that Gq ~ H^Ii as z^/zi+l oo, where M, is the 

limit of Mi{zi, Zn) when Zi/zi-^-i oo. It follows from the theory of difference 
equations that this solution exists and is unique. It can be represented as an infinite 
ordered product of /^-matrices. 

Introduce the matrix-valued function 
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Theorem 4.2 is equivalent to the statement that this function satisfies the difference 
equation 

]C(^zi, Zj, 2;^|s) = 

(4.12) 

Zj-i Z\ Zn Zj+i 



On the other hand, formula (3.24) implies that the function /C satisfies another 
difference equation: 

lC{zi,...,pZi,...,Zn\s) = 
{4:.1S)Gq^{zi, ...,pZi, Zn)Xi{zi, Zn\s)Go{zi, Zn\s)]C{zi, Zi, Zn\s). 

Define the operators Dj acting on ® ■ • • ® as follows: 

(4.14) D,{v, ® • • • ® ^n) = <^^^^j>-<>^->-^^^j-^> ^^^ ^ . . . ^ ^^), 

2{k + 1) 

where Xj are defined by (4.4). From the definition of the fundamental trace it 
follows that J-'{zi, Znis) = J-'o{zi, Zn\s)zi^ ...z^'^, where .Fq — > 1 as s — > 
and Zi/zi^i — > oo. 
Let 

(4.15) e{z\p) = Yl (1 -j9"^^)(l -p'^+'^z-^){l-p'^+^). 

m>0 

Then the function 

is p-periodic: e{pz, a\p) — e{z, a\p). Therefore, the function 

(4.17) K:o{zi,...,Zn\s) = K:{zi, Zn\s)Y[s{z, -Di\p) 

i 

is a single-valued meromorphic function in the region Zi ^ which still satisfies 
equation (4.13). 

Thus, we have proved the following theorem. 

Theorem 4.3. The function /Co satisfies a pair of matrix difference equations - a 
p- difference equation from the left: 

(4.18) K:q{zi, ...,pZj, Zn\s) = Yj{zi, Zn\s)}Co{zi, ...Zj, Zn\s), 

and an s -difference equation from the right: 
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where 

(4.20) Yj{zi, Zn\s) = Gq^{zi, Zn\s)Xj{zi, Zn\s)Go{zi, Zn\s), 



Uj{zi,...,Zn\s) = 

\i ) ^i-l ^1 + 1 

(4.21) ('n£(^,-Ab)). 

We know that the function Uj is a product of f)-periodic functions, so it is p- 
periodic (eUiptic): Uj{zi, ...,pzi, ...,Zn) = Uj{zi, ...,Zi, ...,Zn)- Therefore, we have 

Corollary. The function Yj is s-periodic for all j: 

(4.22) Yj (zi , — , szi, Zfi) = Yj (zi , — , Zi, z^) ■ 

Proof. Since Uj is p-periodic, the function }Co{zi, ...,pzi, Zn) satisfies (4.19) as 

long as /Co(-2i, ^i, ^n) does. Thus, the function /Co(-2i, ■■■,pzi, Zn)JCo{zi, Zi, Zn)~^ 

has to be s-periodic. But this function is exactly Fj, Q.E.D. 

5. Diflference equations with elliptic 
coefficients and double difference systems 

In this section we survey the theory of systems of difference equations, in partic- 
ular those with elliptic coefficients, and introduce the notion of a double difference 

system which naturally arises from this theory. 

Let p and s denote two nonzero complex numbers such that < |s| < 1. If 
/(z) = /(^i, Zn) is a matrix- valued function of n complex variables then let 

^^i/(z) = /(^l, ■■■^PZi, ...,2n), 

(5.1) Sif{z.) = /(Zi, SZi, Zn). 

We will say that a meromorphic function / in C*"' is p-elliptic (respectively, 
s-elliptic) if Pif = f (respectively 5"^/ = /) for all i. 

Let ai{z), 1 < i < n he arbitrary meromorphic functions in C*"^ with values in 
N X N matrices. Consider the system of difference equations 

(5.2) P,/(z) = a,(z)/(z), 
where / is an x A'^-matrix valued function. 

Definition. A system of p- difference equations is called consistent if there exists 
a meromorphic solution f to this system whose determinant is not identically zero. 

The following obvious proposition gives necessary conditions for consistency of 
re o^ 
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Proposition 5.1. If system (5.2) has a nonzero meromorphic solution f such that 
det{f) is not identically zero, then 

(5.3) Pi(ij{^) ■ = -PiOi(z) ■ 

Definition. If a system of difference equations satisfies (5.3), it is called holo- 
nomic. 

Thus, any consistent system must be holonomic. 

Remark. It is obvious that if system (5.2) admits a meromorphic solution / 
whose determinant is not identicaUy zero then aU meromorphic solutions of (5.2) 
form an A'^^-dimensional vector space over the field of p-elliptic functions: all of 
them have the form fg, where is a matrix- valued p-elliptic function. 

In the special case when the coefficients are rational functions the above 
necessary conditions are also sufficient: 

Proposition 5.2. Any holonomic system of p- difference equations with rational 

coefficients is consistent. 

Proof. Wc start with proving a technical lemma. 

Lemma. If ai are rational for all i then there exist positive real numbers r^, 
1 < J < ''T', such that the multiannulus A — {{zi., 2;n)|rj|p| < \zj\ < rj, I < j < n} 
does not contain a singular point of af^ for any 1 < i <n. 

Proof of the Lemma. Since ai are all rational, the singular set for the collec- 
tion of functions {af^, 1 < i < n} is an affine algebraic variety in C"^. By our 
agreement it is not the entire C", therefore it is a subset in some algebraic hy- 
persurface prescribed by the equation P{zi, ...,Zn) — 0, where P is a polynomial: 

— "'mi...m„^l ■ ■ ■ ■ 

Let ri be defined according to the rule: = r, rj-i — e^^ . Denote the corre- 
sponding annulus A defined in the statement of the lemma by Ar- Let T be the 
lexicographically highest nonzero term in the polynomial P (i.e it has the highest 
possible degree of zi, and among the terms with that degree of zi it has the highest 
degree of Z2, and so on). Then it is easy to see that limj-^oo ^ 1 (i-e. the 

highest term dominates all the others). Therefore, if r is sufficiently large, the 
polynomial P cannot vanish on Aj. (because T does not). Thus, all the functions 
a^^ are regular in Ar, Q.E.D. 

Now let us prove our proposition. Pick an annulus A satisfying the condition of 
the Lemma. Now make this annulus into an n-dimensional complex torus (abelian 
variety) as follows: identify two points {zi, Zi, Zn) and {zi, ...,pzi, z^) of A 
whenever \zi\ — r^. Denote the resulting quotient space by (it is nothing else 
but the n-th Cartesian power of the elliptic curve Ep = <C* / {z ~ pz}). System 
of difference equations (5.2) can now be interpreted as a gluing condition for a 
rank N holomorphic vector bundle over , and meromorphic C''^- valued solutions 
of (5.2) can be viewed as meromorphic sections of this bundle. But it is known 
from elementary algebraic geometry that any holomorphic vector bundle over a 
smooth projective variety has nonzero meromorphic sections, and the dimension of 
the space of meromorphic sections over the field of rational functions is equal to 
the rank of the bundle. In our case it implies that (5.2) has N C-'^-valued solutions 
which are linearly independent at a generic point. These solutions can be combined 
into an X A^-matrix solution which will have a nonzero determinant at a generic 

r\ 'cr' T\ 
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However, in general a holonomic system of p-difference equations with meromor- 
phic coefficients does not have to be consistent. A counterexample already exists 
for systems with elliptic coefficients. 

Example. Consider the system of difference equations: 

(5.4) f{pz, w) = f(z, w), fiz,pw) = a(z)f(w), 

where a is a non-constant p-elliptic function in C*. It is clear that this system satis- 
fies identities (5.3). Still, let us show that it does not have a nonzero meromorphic 
solution. 

Suppose that / is a nonzero meromorphic solution of (5.4). Then for all w but 
a countable set of them f{z, w) is a meromorphic function in z. Also, the order of 
f{z, w) as a function of z at any fixed point zq is an essentially constant function of 
w - ii has the same value d{zQ) at all points w except for a countable set of them. 
However, the second equation in (5.4) shows that if zq is a pole of a of order k then 

(5.5) ordz^f{z,p^w) = oidz^ f{z,w) - km, 

so for sufficiently large m the order of f{z,p'^w) at zq will become less than d{zQ) 
- contradiction. 

Remark. Of course, one can find plenty of nonzero solutions of (5.4) with 
essential singularities, e.g. 

Q{w\p) 

f{z,w) = 



Q{wa{z)\p) 

where 6 is defined by (4.15). 

In spite of this example, it is very easy to show the existence of meromorphic 
solutions for a single difference equation. 

Proposition 5.3. Let a{z) be a meromorphic N x N -matrix function in C* such 
that det{a) is not identically zero. Then the difference equation 

(5.6) f{pz) = a{z)f{z) 

has a meromorphic solution whose determinant is not identically zero. 

Proof. Let r be a positive real number with the property: the function a{z) is 
defined on the circles \z\ = r and \z\ = |p|r, and its determinant does not vanish 
anywhere on these circles. Let A = {z & C*\r > \z\ > \p\r} be the annulus squeezed 
between these circles, and let Ep = A/ {z ^ pz ii \z\ = r) be the elliptic curve ob- 
tained by gluing the boundaries of the annulus A to each other. Difference equation 

(5.6) can now be interpreted as a gluing condition for a rank N holomorphic vector 
bundle over Ep. The rest of the argument is as in the proof of Proposition 5.2. 

Now consider system (5.2) in which the coefficients Oj are s-elliptic. Then the 
system has a new property: if /(z) is a solution then Sif{z) is a solution as well. 
If det(/) is not identically zero, this property implies that there exist p-elliptic 
matrix- valued functions 6i(z), 1 < i < n, such that 

(5.7) 5,/(z) = /(z)6,(z). 

We see that there is another system of difi^erence equations with elliptic coeffi- 
cients (now s-elliptic) satisfied by /. Thus, we are naturally lead to introduce a 
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Definition. A double difference system is a system of difference equations of the 
form 

Pif{z) = ai{z)f{z), 

(5.8) Sif{z) = f{z)bi{z), i = l,...,n 

where f is an N x N -matrix valued function, and ai, hi are meromorphic N x A'"- 
matrix valued functions whose determinants are not identically zero. 

Let us now study double difference systems. 

Suppose tliat / is a solution of (5.8) whose determinant is not identically equal 
to zero. Then one must have a consistency condition 

(5.9) ai{z)f{z)Pibj{z) = 5,ai(z)/(z)6,(z), 

since both sides of (5.9) are equal to PiSjf{z), according to (5.8). The simplest 
way to satisfy this condition automatically is to set 

(5.10) SjUi = ai, Pibj = bj, 

which is the same as to say that are s-elliptic and bj are p-elliptic. Besides this, 
we have the usual consistency conditions 

Piaj{z) ■ ai{z) = Pjai{z) ■ aj{z), 

(5.11) Sibj{z) ■ bi{z) = Sjbii^) ■ bji^). 

Definition. A system of the form (5.8) satisfying consistency conditions (5.10) 
and (5.11) is called an elliptic double difference system. 

Elliptic double difference systems are exactly those arising from systems (5.2) 
with s-elliptic coefficients. Our discussion shows that they are the most natural 
examples of double difference systems. 

Note that in an elliptic double diffierence system, the coefficients of the s- 
difference equations play the role of monodromy (or connection) matrices for the 
p-difference equations, and vice versa. 

Of course, (5.10) and (5.11) are only necessary and by no means sufficient condi- 
tions of consistency of system (5.8). This is demonstrated by the following propo- 
sition giving a necessary and sufficient condition of existence of a nondegenerate 
solution for an elliptic double difference system with constant coefficients. 

Definition. We say that the numbers p, s are generic if for m,k E Z p^ = if 
and only if m ~ k = 0. We say thatp, s are strictly generic if they generate a dense 
subgroup in C* . 

Proposition 5.4. Let p, s be strictly generic, and let ai and bi be constant ma- 
trices for all i. Then system (5.8) has a meromorphic solution f with det{f) not 
identically equal to if and only if there exist invertible N x N matrices R, L and 
diagonal matrices Mi, 1 < i <n, with integer entries, such that 
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Proof. 

Sufficiency. The function 

n 

(5.13) f{z) = Ll[zf'^R 

i=l 

is a solution of (5.8) whose determinant is not identically equal to zero. 

Necessity. Assume that (5.8) has a mcromorphic solution / whose determinant 
is not identically zero. Then this solution has to be holomorphic and nondegenerate 
everywhere. Indeed, suppose that zq is a singularity of /. Then (5.8) implies that 
for all ki^nii G Z Hil-^i * S^^)'^o is a singularity of / as well. But since p,s are 
strictly generic, the set of these points is dense in C*" - a contradiction which 
implies the holomorphicity of /. Applying the same argument to /~^, we get the 
nondegeneracy. 

Now we need to use a simple lemma from the classical theory of difference equa- 
tions: 

Lemma. The difference equation f{pz) = af{z) has a holomorphic nondegen- 
erate solution if and only if the matrix a is diagonalizable, and its entries are integer 
powers of p. This solution, if it exists, has the form Lz^ where M is a diagonal 
matrix with integer entries, L, R are invertible matrices, and a = Lp^^ L~^. 

This lemma together with the commutativity condition [a^, Uj] — which follows 
from (5.11), implies that the matrices Oj simultaneously diagonalize in a certain 
basis, and their eigenvalues are integer powers of p. In other words, there exists an 
invertible matrix L and diagonal matrices with integer entries Mi, such that 
Oi = Lp^^L~^, and / = Lz^^'^ . . .z^'^R, from which we get (5.12). 

Proposition 5.5. // p, s are strictly generic then the dimension of the space of 
solutions of any double difference system (5.8) (overC) is less than or equal to N"^ . 

Proof. Let zq be a regular point of the coefficients of (5.8) and all their p, s- 
translates. Such a point obviously exists since the singular set has codimension 
1. Then / must be regular at this point, and its value there determines its value 
anywhere else, since the p, s-translates of zq form a dense set in C*"^. 

Remarks. 1. It is easy to show that Propositions 5.3 and 5.4 are true for 
generic p, s which are not necessarily strictly generic. 

2. The dimension of the space of solutions in Proposition 5.5 can be exactly N"^: 
this happens when and hi are scalar matrices. 

3. It follows from Proposition 5.5 that in the case A'^ = 1 and generic s, if 
a double difference system has a nonzero solution, it is unique up to a constant 
factor. 

One of the first major problems in the theory of elliptic double difference systems 
is the consistency problem: 

Consistency problem I. Classify all sets {oj, 6i, 1 < * < n} for which (5.8) has 
a solution whose determinant is not identically zero. 

Another formulation of this problem is: 

Consistency problem II. For a given set of coefficients {a^, 1 < i < n} find all 
possible sets {6^,1 < ^ < n} for which (5.8) has a solution whose determinant is 
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In general, this problem is very difficult, and it is not clear how to approach it, 
even in the one variable case. However, in the case N = 1 (scalar- valued functions) 
it can be solved completely as described below. 

The idea is to explicitly present a nonzero solution to the p-part of the elliptic 
double difference system, and then find the corresponding functions 6^. Then the 
general form of the functions bi is b* — bi^j^, where his a. 7?-elliptic function. 

For simplicity we assume that n = 1; the case n > 1 is analogous but somewhat 
more technical. Our purpose now is to find a nonzero solution of the equation 
f{pz) = a{z)f{z), where a is an s-elliptic function. 

It is known that any scalar-valued elliptic function can be written as a constant 
times a ratio of products of theta functions. More precisely, for every s-elliptic 
function a{z) of one variable there exist a unique constant C e C, A; e Z, and a 
unique finite set ^ C {2; e C*||s| < \z\ < 1} with a function u : A ^ ^\{0}, such 
that 

(5.14) a{z) = Cz^ Y[ e(^/a;|s)"(^); 

and conversely, the function a given by (5.14) is s-elliptic iff Yl/xeA^^-"^) ~ 

Therefore, it is enough to be able to find meromorphic solutions to the equations 

(5.15) fipz) = Cfiz), fipz) = z^f{z), f{pz) = Q{z/x\s)f{z) 

- then a solution of f{pz) = a{z)f{z) can be obtained as a product of such solutions. 
A meromorphic solution of f{pz) — Cf{z) is given by 

(5.16) /(.) - Q'^l"' 



e{Cz\p) ' 

and 

...7. Wfr ^ Oisz\pMCz\p) 

^ ^ "•'^ ' e{Csz\p)e{z\p) 

A meromorphic solution of f{pz) — z^f{z) is given by 

(5.18) f{z) = {-e{z\p))-\ 

and 

(5.19) /(..)//(.) = (IgPl 
Finally, a meromorphic solution of f{pz) = Q{z/x\s)f{z) is given by 



k 



00 00 
(5.20) f{z) = H (l-pyz/x)-' H {1-pWx/z)Q{z\p)Q{cI>{s)z\p)-\ 
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where (p{s) = Y[m>i(^ ~ ^"^) Euler product, and 

e{z/x\p) e{sz\p)e{(p{s)z\p) 



(5.21) fisz)/f{z) 



(j>{p) Q{z\p)Q{s<j>{s)z\py 



The above formulas allow us to construct a solution of f{pz) = a{z)f{z) for any- 
elliptic function a as a ratio of products of expressions (5. 16), (5. 18), (5. 20), which 
helps us solve the consistency problem. 

Let K{p) be the field of p-elliptic functions, and let K{pY be its multiplicative 
group. Also, let K{p)' be the multiplicative group generated by constants and 
the functions z, Q{z/x\p). Let L{p^s)' be the subgroup of K{p)' consisting of 
functions of the form f{z)/f{sz), f G K{p)' . Let K{p, s)' = K{p)' /L{p^ s)', and let 
K{p,sY^K{pY/{L{p, s)'nK{py). 

Let F : K[s)' — > K{p, s)' be the group homomorphism defined by: 

(5.22) C ^ (5.17), z'' ^ (5.19), Q{z/x\s) ^ (5.21), \s\ < \x\ < 1, 

It is easy to check that F restricts to a homomorphism F : K{s)* K{p,s)*. 
Moreover, it is clear that the kernel of F is L{s,py, so F gives rise to a map 
Fq : K{s,py — > K(p, s)*. The map Fq is obviously a group isomorphism. 

Now we can formulate the necessary and sufficient condition for consistency of 
system (5.8). If a e K{p)*, let 4'ps{a) be the image of a in K{p, s)*. 

Proposition 5.6. The system of difference equations 

(5.23) f{pz) = a{z)f{z), f{sz) = f\z)h{z) 

in which a is s-elliptic and b is p-elliptic, is consistent if and only if Fo{(f)ps{a)) = 

Finally let us describe an explicit formula for solutions of double difference sys- 
tems, which applies in the case when p, s are strictly generic. Consider the system 

(5.24) f{pz) = a{z)f{z), f{sz) = f{z)h{z) 

where a, 6, / are x A-matrix valued functions, and a is s-cUiptic, h is p-elliptic. 
Assume that / is a mcromorphic solution to (5.24), and let zq be a point at which 
all matrices a{p'^z), b{s'^z), nijUEl,, are regular and nondegenerate. Let z &C*. 
Let {rrij}, {uj} be sequences of integers such that limj^oo s'^^ zq = z. Such 
sequences can be constructed easily with the help of continuous fraction expansions. 
Then the following formula is valid: 

Proposition 5.7. 

nfc — 1 

(5.25) f{z) = hm TT a{pizo)f{zo) TT b{sho). 

j—>-oo 

j=mk-l j=0 

Existence of this limit and the fact that it equals f{z) follows directly from 
(5.24). 

Formula (5.25) determines f{z) up to a finite number of unknown parameters - 
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Let us now construct another consistent elliptic double difference system related 
to system (4.18), (4.19). 

Consider an elliptic double difference system of the form 

(5.26) Pi/(z,p,s) = ai{z,p,s)f{z,p,s),Sif{z,p,s) = /(z,p, s)fei(z,p, s), 

where ai,bi, f are N x A'"-matrix valued functions whose determinant is not iden- 
tically zero, and let us assume that bi{z,p,ps) — bi{z,p,s). It is obvious that this 
condition is satisfied for system (4.18), (4.19). Assume that / is a solution to the 
system, and let 

(5.27) c(z, p, s) = /(z, p, ps)f{z, p, s)~^. 
Then we can write two equations for c. First of all, 

PiSifiz,p,ps) = f{z,p,ps)bi{z,p,ps) = c(z,p,s)/(z,p,s)6i(z,p,s); 

(5.28) 

PiSif{z,p,ps) = PiSic{z,p,s)PiSif{z,p,s) = PiSic{z,p,s)ai{z,p,s)f{z,p,s)bi{z,p,s), 
which implies: 

(5.29) PiSic{z,p,s) = c{z,p,s)ai{z,p,s)~^. 
On the other hand, 

(5.30) Pif{z,p,ps) = aj(z,p,ps)c(z,p,s)/(z,p,s) = Pic{z,p, s)ai{z,p, s)f{z,p, s), 
from which we have 

(5.31) ai{z,p,ps)c{z,p, s) = Pic{z,p, s)ai{z,p, s) 
Equatins (5.29) and (5.31) together imply: 

(5.32) Sic{z,p, s) = Siai{z,p,ps)c{z,p, s). 

Thus, we have shown that the function c(z,p, s) satisfies a new elliptic double 
difference system (5.29), (5.32), which involves the functions and does not involve 
bi. 

Now assume that as s — > 0, f{z,p, s) has a finite limit f{z,p, 0) which is known 
and generically nondegenerate. (This property holds for system (4. 18), (4. 19): the 
fundamental trace converges to the highest matrix element of the intertwiner as 
s ^ 0). Then one can write the following formula for /: 

oo 

(5.33) f{z,p,s)^Ylc{z,p,p^s)-'f{z,p,0). 

3=0 

Thus, if we knew an explicit expression for c, we could get a more structured 
expression for / than that coming from (5.25), which would also be free from 
unknown parameters. 

Unfortunately, for system (4.18), (4.19) it is not clear how to compute c explicitly; 
however, the quasiclassical limit of c can be understood. This quasiclassical limit, 
i.e. the first term of the Taylor expansion of c near q=l (where / = /Co) is equal to 
the right hand side coefficient of the "moduli equation" for the fundamental trace 
of the classical affine Lie algebra ([E], Eqn. (3.24)) - the equation characterizing 
the derivative of the fundamental trace with respect to the modular parameter of 
the corresponding elliptic curve. This fact is reassuring, since it gives rise to a hope 
that the method of [E] can somehow be generalized to the quantum case, which 
would allow one to produce some kind of explicit expression for the function c, and 
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Appendix: limiting cases 

Let us briefly describe some interesting limiting cases. These cases correspond 
to some special (limiting) values of the parameters q, s. 

Case 1. The quantum KZ limit: s = 0. In this case, the fundamental trace 
transforms into the highest matrix element of a product of intertwiners, system 
(3.24) becomes the quantum KZ system of Frenkel and Reshetikhin, and relations 
(4.9) become the connection relations for the quantum KZ system. In this limit, 
one gets g-hypergeometric functions and their generalizations. 

Case 2. The elliptic KZ limit: q ^ 1, p = q~'^^^^^\ k is flxed. In this case, the 
fundamental trace becomes the fundamental trace for the classical afflne Lie algebra, 
system (3.24) degenerates to the elliptic r-matrix system involving Belavin's elliptic 
r-matrix, and relations (4.9) degenerate to the monodromy relations for the elliptic 
i2-matrix equations (see [E]). In this limit, one gets transcendental functions of an 
elliptic curve, vector-valued modular forms etc. 

Case 3. The Yangian limit: q = e^, p = k is flxed, Zi = e^^\ e — ^ 0. 

In this case, the degeneration of the fundamental trace should be something like 
the fundamental trace for the (doubled) Yangian of sIn, and the equations (3.24) 
should converge to a trigonometric deformation of the Smirnov's equations [cf [LS]] 
This limit is still unexplored. 

Case 4- The critical limit: p ^ 1. In this case, equations (3.24) transform into 
an elliptic analogue of the Bethe ansatz equations. The Bethe ansatz equations are 
obtained if one combines this limit with the limit s — (cf. [TV]). 

Besides these, there are many other unexplored limiting cases which are easier 
to study than the general case. It is expected that studying these limiting cases, 
one should be able to get interesting information about various classes of special 
functions arising in representation theory of Lie algebras and quantum groups. 
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